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Abstract 

In this paper, we consider a couple of alternative dark energy models using the total equation of 
state of the cosmological fluid, wtot- These models are fit to the recent type-la supernovae data and 
are compared to previously considered models. The first model is based on the hyperbolic tangent 
' and provides a good estimate of the rate of the transition to dark energy domination. The second 

■ model is a cubic spline model. This model demonstrates and quantifies the non-monotonicity in the 

. total equation of state coming from the supernovae observations. At present, the supernovae obser- 
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vations indicate significance to non-monotonically decreasing dark energy. We derive constraints 



D \ on the spline paramters and compare and constrast the results to the Cosmological Constant dark 



energy model. Both the hyperbolic and splines models indicate that a precise physical notion of 
dark enegy is a potentially ever more mysterious quantity? 
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I. INTRODUCTION 



Observations of distant tvpe-Ia supernovae have shed great hght on the evolution of 
universe at late times[l|, 0, H 0, 0, @|- These observations indicate that the universe is 
presently undergoing an accelerated expansion that started about 7 billion years ago. The 
source of this acceleration has been labeled dark energy and the corresponding energy density 
is peculiar in that it appears to require a negative equation of state |8|. Cosmo logists were 
quick to pursue models which reproduce the observed effects of the dark energy. All of 
these models have negative equations of state coming from positive dark energy density and 
negative pressure. 

The first of these models has its historical origins with Einstein. It is the Cosmological 
Constant Cold Dark Matter (LCDM) cosmological model. This model adds an additional 
constant A to the Einstein Equations. The resulting model of dark energy is typically 
characterized by a constant equation of state = Pa/ Pa = Initially, the LCDM 

model fit observational data quite well[l|, [3|. More recent data suggests that while this 
model is not ruled out, it is statistically not as favored as models with smaller constant 
equations of statejl, [9|. Still, others propose to introduce phenomenological models of the 
dark energy equation of state (and consequently the energy density) based on kinematics |lOl| 
or linear Taylor's Series expansions of the dark energy equation of state (wy) [l, [TH, Il2j |. 
The latter model was shown to be divergent at early times and has since fallen out of favor. 
However, this inspired others to introduce alternative evolving dark energy equations of state 
which removed the early time exponentially divergent behavior 13|, The first example of 
these models was originally proposed by [l5| as an alternative to constant equation of state 
models and has an evolving equation of state. In [l6|, they use this model to analyze different 
properties of the supernovae data and implications for the dark energy parameters. For dark 
energy, they find that it is difficult to tightly constrain several dark energy parameters using 
the recent supernovae data. 

In j^, they propose several ansatz models for the hubble parameter which are polynomials 
of the cosmological scale factor. In these models, the dark energy is characterized by an 
equation of state that metamorphosizes from wy = at a redshift 2; > 1 to wy ~ — 1 at 
z = 0. From this analysis, they conclude that the supernovae data favors negative evolving 
equations of state for dark energy. Similar ansatz models for dark energy were proposed 
in [l7[] and are consistent with the results from j^. Still others have proposed modeling 



dark energy equation of state by using cubic splines[l8l]. Here, they constrain dark energy 



in a series models and find that they are consistent with negative equations of state using 
supernovae, gravitational lensing and large scale structure. The results of their analysis put 
tight constraints on dark energy equations of state. 

A very curious set of models are the Sudden Gravitational Transition models, see [Ji 
and references therein. These models suppose that dark energy is the result of a late- 
time phase transition in the universe. The source of the dark energy is the result of a 
resonant effect of a free, ultra-low mass, quantized scalar field coming from Quantum Field 
Theory in Curved Spacetime. With a suitably chosen order parameter rj, this resonance 
causes rj — x^m^ = constant where rj takes the form R^, R^^R'^'^ and R^rupR'^^'^'^- The 
model for which rj_= R^ is called the Vacuum Cold Dark Matter (VCDM) cosmological 



model|20l. |21| . |22| . |23| . |24J |. In all of the Sudden Transition models, the effects of the scalar 
field are negligible until about half the age of the present universe, denoted in terms of 
redshift as Zj and is labeled the redshift of transition. Here x^m^ is the one free parameter 
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of these models with x corresponding to a dimensionless parameter fixed by the theory and 
m is the mass of the scalar field. The constant rj arising from this effect reacts back on the 
universe which results in an accelerated expansion through a sort of gravitational Lenz's 
law. It has been shown that these types of models are tightly constrained by the recent 



cosmological experimental data, see [19|, l25|, l26|, |27j for details. 

The effects of these models are different than the previous ones in that they change the 
nature of gravity itself. Dark energy is not an unknown extra quantity that appears in the 
total energy density ptot but is a manifestation of the changes in gravity induced by the 
effects of the field. However, these changes lead to comparable results to the previous dark 
energy models [l9|, 26, 27 . 



A similar model to the Sudden Gravitational Transition models is the Super- Acceleration 



model[28|, |29|. This model is based on Quantum Field Theory in Curved Spacetime as well. 
A fully renormalized energy density and pressure are computed for a scalar field with a 
quartic self-interaction in a locally De Sitter spacetime. The resulting cosmology leads to a 
dark energy term in the Einstein equations which has a dark energy equation of state which 
relaxes gradually from to -1. 

What we propose here is that we model directly the Hubble parameter's evolution through 
the total equation of state Wtot- This is similar to the method from [sl mentioned above. 
This approach has a novelty that it does not have the present day matter density VLmo (or 
alternatively the present day dark energy density VLxo) as a free parameter and models 
directly the kinematics contained in the supernovae data, see js, 1^ for nice discussions. We 
will take two approaches to the modeling. The first approach is based on a phenomenological 
model that was inspired by the hyperbolic tangent, hence forth labeled the hyperbolic model. 
This model has two parameters which determine the time of transition (a) and the strength 
of transition (/?). This 2 parameter model is fit to supernovae data using a finite linear 
grid and marginal parameter estimates are extracted using the procedure discussed in 27 . 
When suitably tuned by this fit, these parameters cause Wtot to undergo a transition that will 
cause it to deviate from at a redshift z ^ 1. This transition is a monotonically decreasing 
function of time bounded below by -1 or conversely a monotonically increasing function of 
redshift bounded above by 0. 

For simplicity, we have set the amplitude of the transition to unity. This sets the cosmol- 
ogy to asymptote to de Sitter spacetime, i.e. Wtot asymptotes to —1. This is motivated from 
the results presented in [isj, where they show that the supernovae data does not constrain 
the dark energy density for redshifts 2; < 0. If some data in the future becomes available to 
constrain the dark energy at these redshifts, then it is possible to adjust the amplitude and 
it would become an additional parameter. We will show that the transition in this model 
is very reminiscent of the models of dark energy discussed above and in particular behaves 
most like the Sudden Gravitational Models and the dark fiuid model[30|. 

The second approach to modeling Wtot that we will take is a phenomenological cubic 
spline analysis. Our interests are to constrain cubic spline models of Wtot to the ranges of 
redshift associated with dark energy domination between redshifts (0 < z < 1). We choose 
two different uniform densities for the splines, consisting of 3 and 6 points. The spline points 
correspond to the model parameters and will be denoted as for i = 1,2,3 for the 3-point 
model or z = 1, ...,6 for the 6-point model. 

We choose to vary the spline points over a finite uniform linear grid. Performing the fit 
to the supernovae data at each point. From this analysis, we conclude that the supernovae 
data favors not just monotonically decreasing Wtot with one extremal point, but allows 
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(and marginally favors for some data sets) transitions which have several extremal points. 
However, implications for monotonically decreasing Wtot shows that the results of the various 
proposed models above are within the constraints allowed by the spline analysis. 

Both of the cosmologies proposed in this paper have an earlier matter dominated epoch 
followed by a later dark energy dominated epoch (i.e. Wtot = 0). In the earlier matter 
dominated epoch, the energy density which would characterize this later epoch appears as 
(presumably) cold (presumably) dark matter like some of the models considered in [8] . Thus, 
these models have some subset of dark matter which undergoes a transition to dark energy. 



This transition is similar to the dark fluid model presented in [30[. There, he supposes that 
the riddle of dark energy and dark matter are one in the same. Dark Matter begins to 
undergo some transition, like a decay, into dark energy at z ~ 1. 

In section [TTl we will propose the two models of Wtot and fit them to the recent supernovae 
data. One is a toy model based on the hyperbolic tangent and the other is a cubic spline 
model. In section IIIH we will discuss comparison of the models considered in the previous 
section to the LCDM model. Also, comparison of the models to the recent WMAP data 
through the shift parameter will be discussed. 

In section llVt we will summarize our results and make concluding remarks about these 
two models and compare them to the LCDM model. It should be noted that implicit in 
this discussion is the Friedmann-Robertson- Walker (FRW) spacetime invariant line element 
given by 

ds^ = -df + a\t) ( + r\de^ + sin\e)d(f') \ , (1) 
\ 1 — fcr^ / 

where k specifies the curvature of the spatial hypersurfaces in spacetime. In this paper, 
we will assume spatially fiat hypersurfaces, i.e. = 0. Also, we assume that the energy- 
momentum-stress tensor is given by a perfect fluid with total energy density {ptot) and 
pressure {Ptot)- 



II. MODELING DARK ENERGY 



Typically, cosmologists interpret the observed late-time acceleration as arising out of 
some non-standard cosmological energy density. In order for this energy density to achieve 
the observed acceleration, it needs to exert a significant negative pressure at late times in 
the universe's evolution. This means that the ratio of pressure to energy density (presumed 
to be non-negative) is negative, i.e. the dark energy equation of state wy{z) = Py/pv < 0. 
All reasonable fitting dark energy models are in agreement on this point. For the LCDM 
model, wv{z) is equal to —1 for all redshifts. For the VCDM model and the other sudden 
gravitational transition model 1^, 25, 27|, wv{z) is shown to be less than —1 for late-times 
after the transition to dark energy domination. Similarly in 13|, [ij, ll8|, they consider a 
significant number of dark energy models and find wy(0) < is favored. 

With this in mind, we propose to model dark energy with no explicit dependence on fi^o 
by using the total equation of state of the universe Wtot given by 



Wtot 



Ptotjz) 
Ptot{z) 



(2) 



where Ptot{z) and ptot{z) are the total pressure and energy density of the universe respec- 
tively. The benefit of this approach of modeling the effects of dark energy is that at late 
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times any parameters introduced into ([2]) which cause Wtot to become negative correspond 
to dark energy parameters. Previously, one had the parameter f2mo with which to contend 
and when combined with other data would put constraints on fitting to the supernovae ob- 
servations. We will constrain the empirically permissible deviations from a Standard Cold 
Dark Matter model. It should be noted that one could assume this approach and arrive 
at all of the expressions for wv{z) arising out of the various dark energy models. So, this 
approach changes nothing for the previous models but allows for a new angle of examination 
of the supernovae observations which can be compared and contrasted to models that have 
already been examined. 

Now, assuming the invariant line element in ([1]), then the time-time Einstein Equation is 
given by 

H\z) _ 8nG pUz) 

-"O "^-"0 Ptotfl 

where Hq is the present day bubble constant and ptot,o is the present day value of the total 
energy density. We are assuming a spatially flat universe. Thus, the total energy density 
Qq = (SttG)/ {3Hq) Ptotfl = 1- Substituting this result into Q gives 



H\z) ^ ptot{z) 

Hq Ptotfl 



(4) 



Assuming T^^'^ is a perfect fluid, VfjT^'^ = implies that ^i^iM satisfies the following 



expression 

Ptot{z) 



Ptotfl 



Exp 



dz' 

;{1 + Wtotiz')) 







(5) 



1 + z' 

where we have explicitly written the redshift dependence of Wtot- 

From this equation, we can see that modeling wtot is directly modeling the bubble pa- 
rameter, (jl]), which is the method used in [8|. So in this paper, we are presenting alternative 
models of the bubble parameter. 



A. The Hyperbolic Model 

Now, lets suppose a particular phenomenological model of Wtot which will characterize 
the transition to dark energy domination. This model is inspired by the hyperbolic tangent 
function and will be labeled the hyperbolic model. We will assume a prior matter dominated 
epoch foYz^l and above, i.e. Wtot = 0. As discussed above in order to get a dark energy 
model with late-time negative pressure, Wtot must become negative as z goes to (present 
day). Let Wtot be defined by 

^ Exp[-/?(2: - a)] 

'^'"^ Exp[-p{z-a)]+Exp[+f3{z-a)]' 

where /3 and a are labeled as transition parameters. (3 controls the strength of the transition 
to dark energy dominated epoch and a controls the time of transition. 
Substituting into we get 



^ = Exp 

Ptotfl 



dz' ^ Exp[-/3(;z' - a)] 



l + z' \ Exp[-/5(2' - a)] + Exp[-|-;3(z' - a)] / . 



(7) 
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FIG. 1: Plots of the wtot and 2a confidence region for the hyperbolic dark energy model given by 
§il for the supernovae data sets R04, R06, SNLS, R04+SNLS and R06+SNLS data sets. The solid 
curves in each graph correspond to the marginal estimate for a for each data set presented in table 
m The dashed lines correspond to the uncertainty in a. For these graphs, we set /? = 3 which is near 
the center of the ranges given in the table. From the graph, this model starts out asymptotically 
in a matter dominated stage,wtot = 0. As z — > 1, the dark energy becomes significant driving wtot 
negative. This model by construction asymptotes to a LCDM model as z —1. 



Substituting this expression into the time-time Einstein equation (jlj) gives the following 
expression for the bubble parameter for this phenomenological model 



^72- = Exp 

^0 



dz' ^ Exp[-/3(z' - a)] 



/o l + z' y ' Exp[-/3(z' - a)] + Exp[+/5(;z' -a)]' ' ^^"^ 

Now, we fit to the experimental data from R04 and SNLS through the luminosity distance 
following the procedure given in [2^, [27| using ([8]) to give the bubble parameter as a 
function of redsbift. 

Performing this fit, we get the results presented in tabled The first column corresponds to 
the data set(s) used in the fitting. The second column corresponds to the local minimum in 
X^- The last two columns correspond to the marginal parameter estimates and uncertainties 
of a and jS. The resulting fits for each data set are marginally better than for the LCDM 
model. Also, the transition and evolution of the dark energy appears to favor a greater 
negative pressure consistent than that which the LCDM model favors. This is shown for 
each data set in figure [H 

From this figure, the transition to dark energy domination {zj) occurs at z ~ 1 and is 
smooth like the LCDM. From the fit to the supernovae data, this redshift is smaller than 
the corresponding redshift of transition for the LCDM model. However, the hyperbolic 
dark energy increases very rapidly and quickly dominates the energy density in the Einstein 
Equation (jlj). So, while the hyperbolic dark energy model near zj is smooth like the LCDM 



the transition of the dark energy is more akin to the Sudden Gravitational Transition [19 



and Dark Fluid Model 30 



From table [H the a parameter is constrained and (3 appears weakly constrained. However, 
the trend from this table is that the most recent supernovae observations are beginning to 
provide additional information which results in tighter constraints on cosmological parame- 
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TABLE I: A table of the marginal parameter estimates of the hyperbolic model derived from the 
R04, R06 and SNLS data at 2cr. 

Data Set Xmm ^ 

R04 175 0.25 ±0.08 7.0l?:[j 

SNLS 62 0.37 ±0.10 2.1+|^ 

R04+SNLS 238 0.29 ± 0.08 3.4 ± 3.2 

R06 201 0.26 ±0.08 6.3 ± 5.0 

R06+SNLS 264 0.29 ± 0.07 3.2 ± 2.8 



ters. Thus, the supernovae data only constrains one of our transition parameters. In figure 
dl we see a plot of the 2a confidence region for Wtpt coming from the marginal estimates 
and uncertainties. This is in agreement with [l^, 13, 18, [l^, 27|, where they show from the 
analysis of many models that only one parameter of dark energy is constrained by the data. 
This also a gree s with the results obtained below in section IIIII 

In (13I . Il5[|. their analysis assumed particular forms of the dark energy density 
{pv{z)/ pv{0)) with the constrained parameters roughly corresponding to the dark energy 
equation of state wv{z = 0) and w'yiz = 0). One of the simplest model is of the form 

2; 

wv{z) =wo + wi—— , (9) 
i. -\- z 



where Wq and Wi are both parameters of the model[15|. This model has been labeled as 
model 2.0 in [131. We will assume this label here. 



In [18[, they used spline approximations of the dark energy equation of state. We will 
consider a similar model below. In [13] , he assumed a kinematic expansion of the cosmolog- 
ical scale factor up to the snap term (oc a{t)). The Sudden Gravitational Transition models 
from fiol] assumed a spatially fiat universe where the one free model parameter corresponds 
to Qmo- This parameter was found to be tightly constrained. However, if one introduces 
curvature, Q^n is found to have a high degree of covariance with Quo which is not tightly 



constrained 27]. This result is similar to that shown in figure 4 of (igI ]. 

A comparison plot of several of these models and the hyperbolic one is given in figure [2l 
From this figure, we see that all of the models appear to converge in their prediction of Wtot 
at z ~ 0.19. Thus at this redshift, the supernovae predicts that Wtot ~ —0.6. Notice that 
the Linder 2.0 model and the Linear Taylor's Series model have divergences. The former 
diverging in the distant future of z = — 1 and the latter diverging at 2; = — 1 and 2; — > 00. 
The hyperbolic model proposed here has no such divergences and most closely resembles the 
VCDM model and references therein) and the dark fluid model of (30| . 

From the early redshift behavior, it would appear that this model suffers from a similar 
divergence problem as the linear Taylor's Series model used in [s], [13] , but this is not the 
case. Upon inspection of ([5]), one sees that Wtot vanishes for z ^ 1 and greater. This means 
that during the preceding matter dominated stage, this form of dark energy behaves just 
like pressureless non-relativistic matter, i.e. grows as (1 ± zY and has an equation of state 
Wy = 0. Thus, this dark energy model is some form of the dark fluid from S^] and the 



ansatz models from [8]. For z > 1, Qy{z){= 8tt G pv{z) / {3 Hq)) is given by 

nviz)^{Qvj-^mo){l + z)\ (10) 
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FIG. 2: (color online) A comparison plot of Wtot for various dark energy models. The light blue 
curve corresponds to the linear Taylor's Series model from R04, with parameters = 0.27, 
Wo = —1.3 and w'q = 1.5. The green curve corresponds to the hinder 2.0 model (defined by ([9])) 
with parameters wq = —2.25 and wi = 0.3. The purple curve corresponds to the N=4 VCDM 
model with 0^,0 = 0.45. The dashed curve corresponds to the LCDM model with fimo = 0.27. 
The solid black curve is the phenomenological model presented in this paper with a = 0.25 and 
P = 3. 



where flyj is some constant determined by Wtot at the redshift of transition (zj) which is 
a function of f3 and a. For typical values of these two parameters coming from fitting to 
the supernovae data (below), Qvj ~ l-5f2mo- Thus at early time, this form of dark energy 
would only make the universe appear that it has more cold matter than it really does. This 
is the essence of the dark fluid model where dark energy arises out of dark matter around 
the transition at z = Zj. We can compare the hyperbolic model to the other dark energy 
models by introducing the assumption that ptot = Pm + Pv where pm and pv are the energy 
densities of matter and dark energy respectively. Upon inspection of (Q and ([5]), pv for 
the hyperbolic model is present at all stages of the universe's evolution, which is a trait 
that it has in common with the LCDM model. To make a quantitative accessment, it is 
necessary to assume a value for the ratio of the matter energy density to critical density at 
present day, flmo = 87rGpmo/(3-ffo)- Since the hyperbolic model most closely resembles the 
VCDM model, we will assume that Qmo = 0.45 which is the marginal estimate from that 
model. As shown in figure [3], the dark energy density associated with the hyperbolic model 
decreases as z approaches 1 from above behaving as matter, pv oc {1 + z)^. At 2; ~ 1, the 
dark energy evolution begins to deviate from (1 + z)^ and begins to grow in significance. 
This is where the negative pressure from the dark energy begins to grow resulting in the 
observed acceleration. This is indicated by the increase at low redshift in the figure. The 
behavior of this equation of state is reminiscent of the Sudden Gravitational model with a 
non-zero cosmological constant [loj and the dark fluid model [30]. 

In this section, we have shown that the hyperbolic model fits to the supernovae data 
quite well as compared to the other cosmological models. In the next section, we generalize 



the approach of using wtot by using spline approximations, similar to those used in [18 
We have fixed by construction that this model asymptote to a LCDM model as z goes to 
— 1. However, one could introduce a non-unit amplitude in the numerator of ([6]) and have 
it asymptote to any value. However, introducing this parameter would increase this model 
to a 3 parameter model. Given that the supernovae data constrains only one of the two 
parameters, introducing a third would most likely not the statistical significance of the model 
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FIG. 3: A comparison plot of the predicted ratio dark energy density to its present day value 
versus redshift (upper) and the dark energy equation of state (lower) for the phenomenological 
model with a = 0.29 and /? = 3.0 assuming il-mo = 0.45 (solid curves) and the LCDM model 
(dashed curves). At increasing redshift this model indicates that the energy density of the universe 
is behaving identical to that of non-relativistic matter. Approaching present day, this decreasing 
energy density reaches a minimum a begins to increase again. This is reminiscent to the effects 



which occur in the Sudden Gravitational Transition models 19|, the dark fluid model proposed 
in [s^] and the ansatz models from Q]. Clearly, the data allows for physical models which have 
significantly different dynamics. 

to the supernovae data. 



B. Spline Approximation 

In this section, we will use cubic splines to approximate Wtot- This is similar to the 
analysis performed in [3]. There they modeled the dark energy equation of state wv{z) 
using cubic splines. The analysis in this paper will assume a prior matter dominated stage 
ending at 2; ~ 1. Thus, we will assume that the spline points of Wtot for redshifts z > 1 
vanish. The model parameters of these spline models are the spline points that specify 
the spline functions for 2; < 1. Recall from the hyperbolic and other dark energy models 
considered in section III At we found that the dark energy effects become significant near 
the redshift of transition, Zj. We want to use the splines to try and constrain what types 
of transition and evolution of the dark energy from Zj{^ 1) to present day. Thus, the 
redshifts of the parameters of the spline model will be contained in the redshift interval 
[0,1]. In this interval, we will analyze uniform grids of spline points. Consider two densities 
of spline points, 3-point and 6-point, which will be denoted as a, for i=l-3 or 1-6. Our main 
interest here is constraining wtot using the supernovae data, inferring properties of Wtot and 
comparing the results to other models. We will show that the computed confidence regions 
of the splines models will have significant intersection with most of those from the dark 
energy models discussed in section III A[ 

Assuming that Wtot is given by a cubic spline function, we can fit the resulting cosmological 
models using the same procedure that was used for the hyperbolic model. Fitting to the 
data, we find the minimum and marginal estimates for the 3 point spline is present in 
table ini and for the 6 point spline is present in table IIIII Comparing these two tables, we 
find that there is little statistical difference between the significance levels of the marginal 
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TABLE II: Marginal Parameter estimates at 2a and the minimum for the 3 point sphne model 
with parameters oi, 02 and 03. 



Data Set 



01 



02 



03 



R04 
SNLS 
R04+SNLS 

R06 
R06+SNLS 



174 
62 
238 
200 
262 



-0.5 
-0.3 



+0.5 
-0.6 
-1-0.3 
-0.4 



-0.4 ±0.2 
-0.5 ±0.2 
-0.4 ±0.2 



-0.4 ±0.2 
-0.6 ±0.1 
-0.6 ±0.1 
-0.4 ±0.2 
-0.6 ±0.1 



-1.3 ±0.5 
-0.7 ±0.4 
-0.8 ±0.3 
-1.3 ±0.5 
-0.7 ±0.4 



TABLE III: Marginal Parameter estimates at 2cj and the minimum for the 6 point spline model 
with parameters ai, 02, 03, 04, 05 and ag- We find that the either or both supernovae data sets have 
wide marginal distributions for the the parameters a\ and 02 and thus does not offer significant 
constraints. Subsequently, we will not list them in the table below. As well, we find that the R04 
data does not constrain 03. The SNLS data constrains only two parameters 04 and 05. Together, 
both data sets constrain 03, 04, 05 and ag. Unconstrained parameters are denoted by dashes in 
the table. All uncertainties are at 2cj. Notice that the R06±SNLS minimum is greater in the 6 
point model than in the 3 point model. This is due to parameter space grid is more dense than in 
the 3 point model. Computational time was too prohibitive to fill in the grid to the same density 
as was done for the 3 point model. 



Data Set 



03 



04 



as 



06 



R04 
SNLS 
R04±SNLS 

R06 
R06±SNLS 



174 
61 
237 
200 

264 



-0 3+°-^ 
-0 2+°-2 
-0 2+°-2 



-0.3: 
-0.5 



-0.3 
-0.5 
-1-0.5 
-0.6 

-0 4+°-^ 

-0 1+°-^ 
"-'•■'--0.4 

-0 3+°-3 



-0.6: 
-0.6 



-0.4 
-0.4 
-1-0.4 
-0.4 

-0 6+0-'^ 
-0.7 ±0.4 
-0.7 ±0.2 



-0.9 



+0.5 
-0.5 



n n-HO.6 
'-'■^-0.6 

-0 3+°-3 



-0.3 



-0.3 
-0.4 



3 and 6 point cubic splines. For the 3 point splines, all of the parameters are constrained. 
With the 6 point splines, the R04 data constrains 3 parameters, 04, as and a%. The SNLS 
data constrains 2 parameters, 04 and as. For the R04±SNLS data, we find that 03, 04, as and 
Oe are constrained. With the introduction of 21 new supernovae observations at redshifts 
z > 1, we find that the no new additional spline parameters are constrained. However, 
the uncertainties in the parameters appear to have been significantly decreased. Multiply 
constrained spline parameters is in contrast to most of the models considered in section III AI 
where only 1 parameter was constrained for each model. This is not unexpected for the 
spline models since multiple spline points constrain wtot in the regions where the data is 
most sensitive. The fits obtained here are marginally better than that which is obtained for 
the previously discussed models. 

The supernovae data is given in terms of the difference between the apparent magnitude 
(m) and the absolute magnitude (M). A plot of this difference modulo an open vacuum 
(A(m — M)) versus redshift for the 3 and 6 point marginal estimate obtained from the 
R06±SNLS data fit is given in figure H] along with the VCDM model, the LCDM model and 
the hyperbolic model. Going from the 3 point marginal spline to the 6 point marginal spline 
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FIG. 4: (color online) A plot of the Distance Modulus A{m—M) versus redshift z for the R06+SNLS 
best fit wtot for the 3 and 6 point splines. The dashed line corresponds to the LCDM model with 
^mo = 0.27. The green curve corresponds to the Linder 2.0 model with parameters wq = —2.25 and 
wi = 0.3. The blue curve corresponds to the best fit VCDM model from [19I] with = 0.45. The 
6 point spline model appears to be midway between the LCDM and VCDM models. The 3 point 
spline most closely resembles the VCDM model. Prom this result, it appears that the supernovae 
data favors models with a stronger transition to dark energy domination than that predicted with 
the LCDM and other models. 

in this figure, we find that there is convergence to the marginal curve of the VCDM model. 
This implies that as the spline model favors a strong transition to dark energy domination 
occurring at a redshift z ~ 0.8 which is more like the transition coming from the VCDM 
model. This is indicated by the larger A(m — M) in figure H] which corresponds to greater 
observed dimming of the supernovae. Weaker transitions (like dark energy stemming from 
a non-zero cosmological constant) would generate brighter supernovae and those correspond 
to a smaller A(m — M). Thus, the supernovae data appear to favor models which have 
greater transition than a cosmological constant. 

A plot of the marginal estimate of Wtot and 2a confidence limits obtained by fitting to the 
data is given in figure [5l Also, this figure exhibits the non-monotonic functional behavior 
of the splines. Thus, these spline models indicate that the data permits different functional 
behavior than that which was considered in section III A[ Many of the spline models shown 
in figure [5]have equations of state Wtot which experience non-trivial bouncing behavior in the 
dark energy dominated epoch at a redshift of 2; ~ 0.5. Accounting for parameter covariances, 
we find that such bounces in Wtot are within la confidence for each data set. For the LCDM 
[sl, 01, VCDM [l9| or the hyperbolic model, we find that only one parameter is constrained 
(again assuming spatial fiatness) and each of these models have a monotonically decreasing 
equation of state for the dark energy and consequently Wtot- Also, these models assume that 
Wtot asymptotes to a constant value. While the splines do have significant overlapping of 
confidence regions with the other dark energy models from section III At the resulting wtot 
are not constrained to be monotonically decreasing. 

In [13], model 2.0 has a monotonically decreasing equation of state with two free parame- 
ters {wq and Wa) which are found to be constrained. However, as mentioned previously this 
model has an asymptotic divergence in its dark energy equation of state, i.e. it does not 
asymptote to a constant. Thus for any proposed dark energy model, the supernovae data 
can generically can constrain at most 4 degrees of freedom depending on the assumptions 
of the underlying model. For the 6 point spline, constraints for the R06+SNLS data given 



11 



Or- 
-0.25 
-0.5 
-0.75' 
Wtot -1 
-1.25^ 

-1.5 
-1.75/ 



0.6 0.8 1 

redshift z 




0.6 0.8 1 

redshift 







-0.2 
-0.4 
-0 . 6 ~ 
-0 . 8 ► 

-1 
-1.2 
-1.4 



0.2 0.4 0.6 0.8 1 1.2 1.4 

redshift z 




.2 0.3 0.4 0.5 0.6 0.7 

redshift z 



FIG. 5: (color online) A plot of the marginal estimates and 2a confidence limits of wtot for the 3 
and 6 point spfine models. The results are for the R04, SNLS, R04-FSNLS, R06 and R06+SNLS 
data sets. The top two rows corresponds to the 3 point spline and the bottom two rows corresponds 
to the 6 point spline. The splines show a different functional behavior of wtot than that which was 
considered in the models from section III Ai The functions here can have inflection, local maxima 
and minima. From these plots, we find that the non-monotonic behavior in the dark energy occurs 
as high significance level. Also, it is apparent that the most recent data release of R06 has increased 
this level. 



in table IIIII show 4 parameters being constrained. This represents an improvement over 
the R04 data sets. It is expected that these constraints will be reduced further with the 
upcoming ESSENCE data 

As with the hyperbolic model, consider ptot = Pm + Pde and determine the properties of 
the dark energy coming from the spline model that are present in the observational data. 
Substituting this expression for ptot in dSD and solving for p^E gives 



Pde 
Pdeo 



1 



1-n. 



-Exp 



mO 



dz' 
l + z' 



[l+Wtot{z')) 



mO 



1-n. 



:i + z) 



(11) 



mO 



A plot of this expression is shown in figure [6] for ^2^0 = 0.30,0.45 and 0.60 assuming the 3 
and 6 point spline models. In this figure, the dark energy is behaving similarly to the dark 
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FIG. 6: Plots of the ratio of dark energy density {pDE)to present day value {poEo) for the 3 point 
spline (right) and 6 point spline (left) best fit models. For z > 1, both spline models behaves 
as ordinary non-relativistic matter with wtot ~ 0. At about z ~ 1, the dark energy under goes 
a transition which causes it to differ from ordinary non-relativistic matter. This transition is the 
source of the dark energy dominated epoch in both of the spline models. The curves going from the 
top curve to the lowest curve correspond to QmO = 0.30, 0.44, 0.60 respectively. For QmO > 0.44, 
the ratio goes negative. This implies that the dark energy density can violate the Weak Energy 
Condition (WEC), but the total energy density never violates this condition. Notice for the solid 
black curve (r^mo = 0.44) that the energy density vanishes for z ~ 0.8. In the case of the 6 point 
spline, we find that the dark energy is approximately until about a redshift of 0.9 and rises 
significantly at late times. This is dynamically very similar to the Sudden Gravitational Transition 



models 



19|,l25|. 



energy from the hyperbolic model. Here, the dark energy at 2; > 1 is behaving as matter until 
reaching a minimum at 2; ~ 0.8. From there, its evolution begins to deviate from matter 
and this results in a decreasing wtot and gives rise to the effects that are attributed to dark 
energy. This is not unlike what occurs with previously considered models of dark energy. 
From figure O we see that the decomposition of ptot implies a negative pde for Q^no > 0.45. 
This implies that the dark energy density could violate the Weak Energy Condition (WEC), 
However, the total energy density of the universe is always positive and thus 



WEC is not violated. This could lead to subtle effects on the vacuum solution of Einstein 
Equations similar to those discussed in [s^] for f{R) gravities and might be worthy of some 
future study. It is worth pointing out that the critical value Qmo = 0.45 corresponds to the 
marginal estimate of the VCDM model. 

All of the models discussed in section Hi Al had monotonically decreasing WtotS>. There 
are as of yet no models which possess this bouncing characteristic on scales of order unity in 
redshift. In [24], they show that the VCDM model possesses a bouncing characteristic with 
respect to the order parameter around the time of transition to dark energy domination. 
However, they found that this behavior experiences a rather rapid exponential decay in the 
dark energy epoch and asymptotes rather rapidly to a constant value on a time scale much 
smaller than the bubble time. 

Despite the dynamical differences in wtot between the dark energy models consider in this 
and the previous sections, we show in figure [7] that the various cosmological models have 
significant overlap in or are very close to the 2a limits obtained from the spline analysis for 
both 3 and 6 point splines. For the SNLS 3 point spline fit, the very tight constraints on 
a2 result in the greatest difference of the regions. However, the 6 point spline more than 
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FIG. 7: (color online) A plot of the 2a confidence limits for wtot obtained from the 3 and 6 point 
cubic spline models (top and lower plots respectively). The solid and dashed black curves are the 
best fit and 2cr confidence limits for the hyperbolic model. The blue curves correspond to the 
VCDM model best fit (solid curve) and the 2cr confidence limits (dashed curves). For this plot, if 
no constraint was obtained from the data fit then all values between and -1.4 are possible. For 
sake of easy of plotting, unconstrained parameters were assigned the value of 0. Notice, that the 
regions of greatest intersection of the models correspond to precisely the region shown in figure [2] 
around redshifts of z « 0.3 for the 3 point splines and z « 0.8 for the 6 point splines. 

encompasses all of the confidence regions of the other models. This is more than indicative 
of the potential constraints on cosmological models. As we know from the Xmm's from tables 
[IT] and [Tnl that there is little significance difference between the 3 and 6 point spline models. 
Thus, if one is interested in model independent constraints for dark energy arising from 
these two spline models, then the constraints derived from the 6 point spline are a more fair 
estimate. 



III. DISCUSSION 

Now lets consider the implications of the two models considered above and compare them 
with the recent supernovae results of R06. There they published 21 new supernovae at high 
redshifts z > 1.0. The fitting to the data utilizes the Monte Carlo Markov Chain algorithm 
from Lewis and Briddle[38|. They are able to draw significant constraints on the values of the 
present day dark energy density assuming it arises from a non-zero cosmological constant. 
Included in their analysis are potential non-monotonically decreasing dark energy equations 
of state. They say that these models are ruled out due to Occam's Razor stemming from 
the greater number of model parameters. 

We aim to show that the implications of this analysis are not as clear as they may first 
appear. There are certain assumptions on the total equation of state of the cosmological 
fluid which one must assume in order to have a LCDM model in the context of Wtot- The 
resulting implications for a relative model comparison make direct comparisons quite difficult 
to analyze quantitatively. However with a few simplifying assumptions, a model comparison 
that includes all of the models priors can be made. In the following analysis, we want to 
calculate the relative probability of the LCDM versus a 3 point spline model for the R06 
data. We will compute the Occam Factors and relative likelihoods associated with each 



model given the recent supernovae data, see j39| for a good example. It will be shown that 



the supernovae data permits us to constrain at most 3 or at most degrees of freedom of any 
perspective cosmological model. In the context of this paper, we will take the perspective of 
the total equation of state. This is important since we are trying to constrain the potential 
behavior of dark energy and determine the degrees of freedom associated with dark energy. 
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Understanding how to compare relative probabilities of the two models begins with Bayes' 
Theorem and multiplication property of probabilities, i.e. respectively 



P{X\Y,I) = ^^P{Y\X,I) (12) 

P{X,Y\I) = P{X\Y,I)P{Y\I) , (13) 

where P stands for the probability, Y and X are two propositions, I is a universal set. For 
an arbitrary cosmological model, X stands for the model and Y represents the superovae 
data (SNLS and/or R06). 

For the LCDM model, we want to compute P{LCDM\SN, I). This model has certain 
assumptions or propositions that are built in which distinguish it from the Spline models. 
We will analyzed these in terms of Wtot- The first proposition is monotonicity of Wtot, i-e. 
Wtotizi) < Wtot{z2) for zi < Z2- This proposition will be labeled hence forth as A and the 
apriori probability function is P{A\l). Secondly, Wtot asymptotes to a constant (—1 for 
the LCDM model). This proposition will be labeled B and its prior probability function is 
given by P{B\I). Any calculations of relative model probabilities must factor in the apriori 
probability of these propositions when one compares the LCDM model to a model without 
these assumptions. We will find that the associated prior probability functions associated 
with these propositions will be included in the Occam Factor. 

From ([iId, the function P{LCDM\SN,I) can be written 

P{LCDM\SN, I) = ^^^^^PiSN\LCDM, I). (14) 

Now, we can rewrite the last term on the right hand of this equation in terms of propositions 
A and B as follows 

P{SN\LCDM, I) = P{SN\A, B, I). (15) 
Rewriting P{LC DM\I) using (|T3l) and substituting this and ( |T5l) into ( |T4l) gives 

P{LCDM\SN, I) = ^^§I^^^PiSN\A, B, I) , (16) 

where the last term on the right hand side is proportional to the likelihood function that 
is used to compare predictions of the LCDM model to the supernovae data. Also, we have 
assumed that P{A\r) and P{B\I) are not apriori correlated in anyway. 

We are not quite at the LCDM model. Refinement of the last term on the right of (fT6|) 
is required. Proposition B states that the model being considered asymptotes to a constant. 
For the LCDM model, this is assumed to be — 1. In |9|], they considers deviations from this 
value and even makes it a free model parameter. The only difference between these two 
asymptotic models are the prior distributions functions for the constant, labeled c hereafter. 
Computing the marginal likelihood over all values of c, (fT6l) becomes 



P{SN\B,A,I) = I P{SN\c,B,A,I)P{c\B,I)dc, (17) 
Jd 

where D is the domain of c and P{c\B, I) is the prior probability distribution function for c. 
The LCDM model assumes that this function is given by P{c\B, I) oc 6{c +1). In [1^], this 
assumption is removed. 
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Here, we are concerned solely with the LCDM model and its marginal probability given 
the supernovae data. With the prior on c, this leaves the LCDM model with one free 
parameter, A. So, (fT7|) can be written in terms of this parameter and gives 

P{SN\B,A,I) = [ P{SN\A,c= -l,B,A,I)P{A\c= -l,A,B,SN,I)dA , (18) 
Jd 

where P{A\c = —1, A, B, SN, I) is the prior probability distribution function for the free 
model parameter A. So, combining fll8l) and fll6l) gives 

P{LCDM\SN,I) = ^^1§W1 

(19) 

P{SN\A, c=-l,B,A, I)P{A\c = -1, A, B, SN, I)dA . 

This expression is the probability of the LCDM model given the supernovae data alone. If 
we want to compare the LCDM model to any other model. This is the term that we need 
to use. Now, lets do the analogous calculation for the spline model. 

Lets assume for the sake of simplicity the 3 point spline. We want to compute the 
probability of the model given the supernovae data {P{Sp\SN, I)) as we did above for 
the LCDM model. To do this, we follow a similar argument as before. Using (fT2l) . this 
probability can be written 

P{Sp\SN, I) = ^lljj^fiSN\Sp, I) . (20) 

We can express P{SN\Sp, I) as an integral over the parameter space 

P{SN\Sp,I)= / P{SN, ai, a2, a^lSp, I)daida2das , (21) 
Jv 

where Oi, 02 and 03 are the model parameters and V is the volume of the 3-d parameter 
space. Using (fT3l) . this can be written as 

P{SN\Sp,I) = 

, (22) 

Jy P{SN\ai, 02, as, Sp, I)P{ai\Sp, I)P{a2\ai, Sp, I)P{as\ai, 02, Sp, I)daida2da^ 

where P{SN\ ai,a2,a3, Sp, I) is the likelihood function; P{ai\Sp,I), P{a2\ai, Sp, I), 
P{a3\ai, a2, Sp, I) are the parameter prior probability functions. There is no reason prior to 
fitting to the supernovae data to assume that the parameters are correlated. So, assuming 
apriori that the parameters are uncorrelated gives 



P{SN\Sp,I) = / P{SN\ ai, 02, 03, Sp, I)P{ai\Sp, I)P{a2\Sp, I)P{a-i\Sp, I)daida2da^ . 
Jv 

(23) 

For sake of simplicity, lets assume uniform priors for the prior probability functions in (l23ll . 
This equation can be written 



P{SN\Sp, I) = P{ai\Sp, I)P{a2\Sp, I)P{a3\Sp, I) j P{SN\ai,a2,a3, Sp, I)daida2das . 

(24) 
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If we assume that the prior probabihty distribution function P{A\c = —1, A, B, SN, I) 
in (fT9|) is a uniform distribution as well, then we can pull it out front of the integral. From 
(fT9|) and fl24l) . the relative probability of the LCDM model to the Spline model given the 
supernovae data is 

P{LCDM\SN,I) ^ J^P{SN\A,c=-l,B,A,I)dA 
P{SN\Sp, I) Jy P{SN\ai, a2, as, Sp, I)daida2da3 ' 



where OF stands for the Occam Factor [39| given by the ratio of the priors out front on the 
integrals and has the form 

^ PiA\I)PiB\I)PiA\c = -1, A, B, SN, I) 

Pia,\Sp,I)Pia2\Sp,I)Pias\Sp,I) ' ^ ^ 

When fitting to the data, we approximate the posterior distribution function for a general 
cosmological model given the supernovae data as 

P{SN\a,Mode\,I)da = Exp[-x\a)/2] , (27) 



where a represents the model parameters and defined in [6|. 

The ratio integrals in fl25p are the relative marginal likelihood of the LCDM and Spline 
models that is determined by fitting to the data. If the relative probability only involved 
the ratio of the integral terms, then obviously models with more parameters would be 
favored over simpler models. The point of the Occam Factor is to counteract this ratio 
and keep over complex models from being favored simply due to their greater number of 
degrees of freedom. This is the idea considered in [gI] when comparing more complex spline 
models (with 3 or more parameters) relative to the LCDM model. However, this is not the 
only consideration when determining the relative likelihood. The LCDM model has prior 
probabilities in the Occam Factor that involve the assumptions A and B above. It is difficult 
for us to quantitatively determine precisely what these probabilities are. 

As a simple model, consider that these priors are unity and assume that both models 
will initially begin as Standard Cold Dark Matter models (i.e. no dark energy). We expect 
from the WMAP and other cosmological data that Qa ~ 3/4. Since, we are assuming 
uniform priors for the parameters of both models then it reasonable to expect from the tight 
WMAP that the limits off^A is bounded by the interval [0,0.9]. For the spline models, we 
assumed uniform priors Oi, 02 and as from the interval [0,-1.4]. Plugging this into the (126|) 
gives OF ^ 3. The majority of the probability associated with each marginal distribution 
function is located near the minimum in x^- So, we can approximate the ratio of the marginal 
likelihood distribution functions as 

^Exp[-x'(A)/2]rfA ^ Exp[-x2(Ao)/2] ^^g) 



JyExp[-x'^ia,i,a2,as)/2]daida2das Exp[-x^(aio, 020, a3o)/2] ' 

where the subscript denotes best fit value, and we have set the ratio of the volumes of 
parameter spaces to 1. We found for the R06 data that the spline model has a which 
is about 4 lower than the LCDM model. Thus, (|28ll is approximately 0.2. So despite the 
increase number of model parameters associated with the 3 point spline models, the data still 
appears to favor dark energy associated with splines over the LCDM model. Furthermore, we 
can interpret the implications of the resulting fit and conclude that there is some significance 
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to the idea that assuming propositions A and/or B is too restrictive. In this simple model, 
we have OF ^ 3 and the ratio of the marginal likelihoods is given by e~^ °. Now lets turn 
the argument above around and compute the probability of the A and B being true given 
the supernovae data. So, solving (1231) and for P{A\I)P{B\I) and substituting in the 
values from the previous paragraph gives 

-2.0 /i 43\ 

P{A\I)P{B\I) ^ l^-^j ^ 0.37 . (29) 

Thus in the approximation of uniform priors for the LCDM model, we find that there is 
significance to the idea that the physics is more dynamical than what occurs in the LCDM 
model. This does not say that the LCDM model is to disregarded, but it does say that 
at present the SN data can not distinguish between the spline models (multiple parameters 
and non-monotonic behavior) and the LCDM model (one parameter and monotonicity) . In 
fact, the spline models can serve as a test of some of the underlying assumptions of the 
LCDM model. This leaves open the door to some of the more exotic models of dark energy 
discussed in section |T] which are dynamically very different than the LCDM model. One 
could propose to increase the number of parameters, e.g. introducing spatial curvature to 
improve the fit. In 1^, 27 1, non-flat models are shown to improve the fit but only marginally. 
For the LCDM model, this does lead to a significant expansion of the likelihood contours 
due to the high correlation between A and the curvature parameter, Q^o- However, this 
does not really affect the results discussed above and does not account for the WMAP data 
which places very tight constraints on curvature [l^] (similar correlations were found for the 
curved VCDM models in 0,|23). 



Lets now compare the predictions of the Wtot models with the predictions of WMAP[32|. 
To do this, we will use the results of Wang and Mukherjee(4^. They derive model indepen- 
dent constraints on the shift parameter R = y/Qmo Jq'^'"' dz/E{z) where E{z) = H{z)/Hq, 
flmo is the present day matter density and Zrec is the redshift of recombination 1100). 
The models that we are considering in this paper do not have a clearly defined value of Qmo- 
However, one point that both models have in common and would be indistinguishable would 
be for z ^ zj where zj is the redshift of transition to dark energy domination. In this limit, 
both models are dominated by non-relativistic matter. For the LCDM model, this means 
that H^{z)/Hl ^nmoil + zf. 

For the spline models, z ^ Zj means that Wtot ~ 0. Thus, ([3]) and (j5]) can be written 

Comparing this result to the corresponding result for the LCDM model, we can define 
a dimensionless matter density which evolves in a similar fashion in the spline model as 
follows 

n^o = (1 + ^i)"'Exp (3 J Y^^^ + wtotiz))] . (31) 
With this definition, ([3]) can be written as 



p(z)2 ^ J ^mo(l + Zf for Z> Zj I . 

^ ' ^ fi™o(l + %)'Exp [-3 ^(1 + wtotiz))] ioTz<z, ^ ■ ^ ^ 
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Using the results from fitting to the supernovae,we find that Qmo = 0.3 ± 0.05 at 2a. 
This range has significant overlap with the LCDM estimate of Qmo = 0.27 ± 0.02[32]. Using 
this results and the above expression for the shift parameter, we get i? = 1.8 ±0.1 at 2a. In 
13, they find that R = 1.70 ± 0.03 at la. 

In this section, we have compared model implications coming from the new supernovae 
data to that of the 3 point spline model considered in this paper. We have shown that one 
must be careful in the assumptions of any parameterization of dark energy. Also, we have 
shown that the most recent WMAP CMB data can be consistent with models that have 
non-monotonically decreasing total equations of state. 



IV. CONCLUSION 

We have shown in this paper, that the supernovae data offers very tight constraints on 
Wtot- This analysis assumed two very different cosmological models. The first is a 2 param- 
eter phenomenological model based on the hyperbolic tangent. The supernovae data gives 
parameters estimates of a = 0.29 ± 0.08 and (3 = 3.2 ± 2.8 at 2a for the R06+SNLS data. 
As with many other dark energy models, we find that only the parameter a (corresponding 
to the time of transition) is significantly constrained leaving jS relatively unconstrained. By 
varying these parameters, we find that this model is able to reproduce similar transitions 
to dark energy domination as several previously considered dark energy models. Of these 
models, the hyperbolic model best reproduces transitions like those coming from the Sud- 
den Gravitational Transition Models [l9| and is more physically attuned to the dark fluid 



model 30 



The other model of dark energy that we considered was a cubic spline model of Wtot- 
This model supposes a prior matter dominated epoch for z ^ 1 and has the potential of 
some non-trivial dark energy behavior for z < 1. Two different uniform density of spline 
points were considered one with 3 spline points and 6 points from the interval z G [0,1]. 



These models are similar to those considered in [18j but there they were used to constrain 
the dark energy density. In this paper, we model the total equation of state Wtot and find 
that the supernovae data has a very wide flexibility in the total equation of state. As many 
as 4 spline points can be constrained by the supernovae data with a minimum number of 
assumptions. There is significant overlap of the 2a confidence regions between these models 
and the monotonically decreasing models considered in section III A[ This is consistent with 



results presented in [18|. Also, the models permit the possibility that the dark energy by 
itself can violate the Weak Energy Condition (WEC), but together with matter there is no 
violation of this condition. 

We find that there is general agreement between the cosmological models considered in 
this paper and the cosmological data coming from supernovae and WMAP. At present, 
there is not enough information to determine the precise physical behavior of dark energy 
and thus we have little clue to its origins and/or underlying physics. This work suggests 
that our scope of investigation into dark energy should not be limited to the LCDM model 
or models which are solely monotonically decreasing. The present data permits and even 
marginally favors models with non-trivial dynamics. Perhaps, future observations with 
more accurate observations would give increased knowledge about the dynamical behavior 
of Wtot- We have shown that the update of R04 with the R06 data has reduced the estimated 
parameters uncerainties by about 10%. The upcoming ESSENCE project [il] results may 
provide some additional upper limits on this type of behavior. The proposed Supernova 
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Acceleration Probe (SNAP)j4J] at roughly 2000 observed supernovae per year should be 
able to distinguish between monotonic and non-monotonic Wtot through its observation of 
the time evolution of the equation of state. 

Future work could include determining the precise solutions of the perturbation equations 
which could potential constrain Wtot- This would allow comparison with the growth of Large 
Scale Structure. Examining the subtle effects of these equations of state to post-Newtonian 
approximation may lead to some constraints as well. 
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